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1 Introduction 



This paper, a continuation of |16|, connects the analysis of the length of the longest weakly 
increasing subsequence of inhomogeneous random words to a Riemann-Hilbert problem and 
an associated system of integrable PDEs. That such a connection exists is not so surprising 
given the fundamental work of Baik, Deift and Johansson connecting the related problem 
involving random permutations to a Riemann-Hilbert problem. For the reader's convenience 
we first summarize some of the results of [|l6[ before presenting our new results. 

A word is a string of symbols, called letters, which belong to an ordered alphabet A of 
fixed size k. The set of all such words of length N, >V(.4, A^), forms the sample space in our 
statistical analysis. We equip the space W(^, N) with a natural inhomogeneous measure by 
assigning to each letter i £ A a probability pi and defining the probability measure on words 
by the product measure. We also order the pi so that 

Pi>P2>--->Pk 

and decompose our alphabet A into subsets Ai, A2, . . . , .Aa/, M < A;, such that pi = pj if and 
only if i and j belong to the same Aa- 
Let 

w = a\a2 • • • ctN £ VV, Oi S .4, 
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be a word. A weakly increasing subsequence of the word w is a subsequence ajjOij • • • ctj^ such 
that ii < i2 < • • • < im and < < • • • < The positive integer m is cahed the length 
of this weakly increasing subsequence. For each word w W we define In{w) to equal the 
length of the longest weakly increasing subsequence in The function 

In ■■ W{A,N) ^ R 

is the principal random variable in our analysis, and the corresponding distribution function, 

F7v(n) := Prob {£n{w) < n) , 

is our principal object. (Prob is the inhomogeneous measure on random words; it depends 
upon and the probabilities pi.) 

To formulate the basic result of define 

k — \ A \ 

where 

M 

A — Aa 

is the decomposition of the alphabet A introduced above, then 

(27r)-('^-i)/2 -q^^l; 2! . . . (yt^ _ x 

a 

/ • • • / n Aa(6' e~^^'/^ 6iJ2 e.) dCi • • • d^k, 

where Aq(^) is the Vandermonde determinant of those with i G Aa, and H denotes the set 
of those that .^i+i < whenever i and i + 1 belong to the same Aa- 

This result has the following random matrix interpretation. The limiting distribution func- 
tion (as N oo) for the appropriately centered and normalized random variable ijsf is related 
to the distribution function for the eigenvalues in the direct sum of mutually independent 
ka X ka Gaussian unitary ensembles]^ conditional on the eigenvalues satisfying J2 \/PiCi = 0- 
In the case when one letter occurs with greater probability than the others, this result implies 
that the limiting distribution of {In — Npi)/^/ N is Gaussian with variance equal to pi(l — pi). 
In the case when all the probabilities are distinct, we proved the refined asymptotic result 

E(^^) = Npi + Y, --^ + O(^), ^ oo. 

Pi - Pj VN 

^There may be many subsequences of w that have the identical length £n{w). 
basic reference for random matrices is Mehta's book [ESf. 



lim Prob ( ^ , < = 
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The derivation of the above asymptotic formulae follows from a direct asymptotic analysis 
of the right hand side of the basic combinatorial equation, 



Fioh{£N{w)<n) = ^ sx{p)f 



Here X ^ N denotes a partition of A^, sx{p) is the Schur function of shape A evaluated at 
P ■= {p1jP2, ■■■,Pk,0,0, ...), and equals the number of standard Young tableaux of shape A, 
see, e.g. [|^. After ||l^ was written, Stanley [^] showed that the measure Prob({A}) := sa(p)/'*' 
also underlies the analysis of certain (generalized) riffle shuffles of Bayer and Diaconis 
Stanley relates this measure to quasisymmetric functions and does not require that p have 
finite support. (Many of our results generalize to the case when p does not have finite support, 
but we do not consider this here.) The measure considered here and in is a specialization 



of the Schur measure Prob({A}) := sx{x)s\{y) For the Schur measure, Okounkov [25| has 
shown that the associated correlation functions satisfy an infinite hierarchy of PDEs; namely, 
the Toda lattice hierarchy of Ueno and Takasaki [^2[. Similar results were also obtained by 
Adler and van Moerbeke [01, 531. 



Gessel's theorem |13] (see also [|l6|, ^) implies that the (exponential) generating function 



of Pioh{£]\[ < n) is a Toeplitz determinant^ 

oo ,N 

Gi{n;{p,},t) := ^ Prob (£jv(^x;) < ^) Try = Dnifi), (l-l) 

N=0 



where 



fi{z) = e'/^ 1[{1+P, 



z . 



Probabilistically, Gi{n; {pi},t) is the Poissonization of ^tv- Similar Poissonizations have proved 
crucial in the analysis of the length of the longest increasing subsequences in random permu- 
tations 1^, 1^, § (see also ||2l|, |3^ and references therein). 

In the present paper we use ( |1.1[) to express Gj{n; {pi},t) in terms of the solution of a certain 
integrable system of nonlinear PDEs. Indeed, we show that Gj{n;{pi},t) can be identified 
as the Jimbo-Miwa-Ueno T-function [18, corresponding to the (generalized) Schlesinger 
isomonodromy deformation equations of the 2x2 matrix linear ODE which has M + 1 simple 
poles in the finite complex plane and one Poincare index 1 irregular singular point at infinity. 
Recall that the number M is the total number of the subsets Aa C A. The poles are located at 
and —Pi^ {id = max Aq)- The integers ka appear as the formal monodromy exponents at the 
respective points —pi^ ■ We also evaluate the remaining monodromy data and formulate a 2 x 2 
matrix Riemann-Hilbert problem which provides yet another analytic representation for the 
function Gj{n; {pi},t). Similar to the problems considered in and j^, the Riemann-Hilbert 

^If (f> is a function on the unit circle with Fourier coefficients 4>j 1/2it J^^ e~^^^ 4>{^^^)dd then T„(4>) denotes 
the Toephtz matrix {(pi-j)^ ^ and Dn{(j)) its determinant. 
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representation of G j {n; {pi} , t) can be used for the further asymptotic analysis of the random 
variable £n{w) via the Deift-Zhou method In the homogeneous case, i.e. when M = 1, 
the system of Schlesinger equations we obtain reduces to a special case of Painleve V equation. 
This result was obtained earlier in . The exact formulation of the results indicated above 
is presented in Theorem 1 in §4. 

Our derivation of the differential equations for the function G j {n; {pi} , t) follows a scheme 
well known in soliton theory (see e.g. |2^]) called the Zakharov-Shabat dressing method. We 
are able to apply this scheme since there exists a matrix Riemann-Hilbert problem associated 



to any Toeplitz determinant as was shown by Deift |10|. For the reader's convenience this is 
derived in §2. 



The basic idea of the Riemann-Hilbert approach to Toeplitz determinants suggested in [10| 
is a representation of a Toeplitz determinant Dn{(j)) as a Fredholm determinant of an integral 
operator acting on L2(C), C=unit circle, and belonging to a special integrable class which 



admits a Riemann-Hilbert representation |15|. Borodin and Okounkov (see also Q for 



a simplified derivation and |22], |^ for a particular case of (p) found a different Fredholm 
determinant representation for Dn{(p). The Fredholm operator in this representation acts on 
l2{{n,n + 1, ...}) which makes the representation quite suitable for the analysis of the large n 
asymptotics of Dn{(t>) (see [21| [22|, Borodin Q subsequently observed that the discrete 
Fredholm representation of Q involves a discrete analog of the integrable kernels and can be 
supplemented by a discrete analog of the Riemann-Hilbert problem. (This is similar to the 



pure soliton constructions in the theory of integrable PDEs |24|.) 

We conclude this introduction by noting that our derivation of integrable PDEs for the 
Toeplitz determinant Dn{fi) can be applied to any Toeplitz determinant whose symbol (j) 
satisfies the condition, 

— log (j)(z) = rational function of z. 

dz 

This is one place where the finite support of p is crucial. It is an interesting open problem, 
particularly in light of |2q], to remove this restriction. 



2 Fredholm Determinant Representation of the Toephtz De- 
terminant and the Riemann-Hilbert Problem 

Let (p{z) be a continuous function on the unit circle C = {|z| = 1} oriented in the counter- 
clockwise direction. Let n G N and denote by Kn{(j)) the integral operator acting on L2{G) 
with kernel 

It was shown in [^] that 

Dn{4>) = det(l - i^„(0)), (2.2) 
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where the determinant on the right is a Fredholm determinant taken in L2{C). (Note that 
Kn{z,z') has no singularity at z = z'.) Equation ( |2.2| ) fohows from the "geometric sum form" 
of the kernel K^, 



n-l 



Kn{z,z') = Y,z' 



fc=0 



2m 



'-^{z')-'-\ 



which shows that the Toeplitz matrix Tn{(f)) is essentially the matrix representation of the 
operator 1 — Kn{4>) in the basis {z'^}_oo<A;<oo- (For more details see p^ .) 

The integral operator Kn{4>) belongs to the class of integrahle Fredholm operators [15, 31, 10|, 
i.e., its kernel is of the form 



KJz, z) 



z — z' 



where 
and 



/(^) = (/l,/2f = (^Mf 



g{z) = {g^,g2Y =(^"",-1) 



2m 



We require, so that there is no singularity on the diagonal of the kernel, 

f{z)g{z)=d. 



(2.3) 
(2.4) 

(2.5) 



An important property of these operators is that the resolvent Rn = (1 — K) ^ — 1 also belongs 
to the same class (see again [15|, [^], |l^ ). Precisely, 



Rn{z,z') 



F^{z)G{z') 



z — z' 



(2.6) 



where 



Fj = {l-Kn)-'fj, Gj = il-K^)-'gj, J = 1,2. 



The vector functions F and G can be in turn computed in terms of a certain matrix Riemann- 
Hilbert problem |15|. Indeed, let us define (cf. [15, the 2x2 matrix valued function 



Y{z)=I - I F{z')gT{z')- 



dz' 



z4C. 



c 



z' — z 



(2.7) 



Let y±(^z) denote the boundary values of the function Y(^z) on the contour C, 

y±(z) = lim y+(z')- 

z'siside 



From (PTD it follows that 



y+(z)-y_(z) =-27riF(z)/(z) 



(2.8) 
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and hence (recall (|2.5D) 

Y+{z)f{z)=Y^{z)f{z). 

Using this, the matrix identity 

F{z')g^ {z')f{z) = f'^ {z)g{z')F{z'), (associativity of the matrix product) 



and (2.7), we have 



Y^{z)f{z)=f{z)- I f^^)g^^')F{z')4^ = f{z)+ I K{z,z')F{z')dz'. 
JC z' — z Jc 

From the definition of F it follows that 

F{z)=Y±{z)f{z). (2.9) 
This and (|2.8|) imply the jump equation, 

Y_{z) = Y+{z){I + 2Tiif{z)g{z)), zeC. (2.10) 

This equation, supplemented by the obvious analytic properties of the Cauchy integral in (^.7|), 
shows that the function Y[z) solves the following 2x2 matrix Riemann-Hilbert problem: 

• Y{z) is holomorphic for all z ^ C, 

• y(oo) = /, 

• Y_{z) = Y+{z)H{z), zeC, 
where the jump matrix H is 

H{z) = I + 27rif{z)g{z) (2.11) 

(2.12) 

These analytic properties determine Y uniquely. To see this, we first observe that det H(z) = 1 
implies that the scalar function deiY(z) has no jump on C; and hence, it is holomorphic and 
bounded on the whole complex plane. This together with the normalization condition at z = oo 
implies that det Y{z) = 1. Suppose that Y(z) is another solution. Since both the functions Y{z) 
and Y{z) satisfy the same jump condition across the contour C, the matrix ratio Y{z)Y^^{z) 
has no jump across C. This means that Y(z)Y^^{z) = constant, and from the condition at 
z = oo we actually have that Y{z)Y~^(z) = L The uniqueness now follows. 

Since Y is the unique solution of the Riemann-Hilbert problem, one can now reconstruct 
the resolvent R using (|2.9| ) and the similarly derived identity 

G{z) = {Ylr\z)9(z), (2.13) 
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for G. We shall refer to this Riemann-Hilbert problem as the y-RH problem. 

Following jl^ the y-RH problem can be transformed to an equivalent Riemann-Hilbert 
problem which is directly connected with the polynomials on the circle C orthogonal with 
respect to the (generally complex) weight (j){e^^). To this end we first note that since the entries 
of / are polynomials in z, ( |2.9| ) implies that F is an entire function of z. Since Y{z) I as 
z — > oo, it follows in fact that F is polynomial, 

F{Z) = ^ , Pn{z) = Z" + Qn-liz) = Qn-lZ^'' + (2.14) 

for some constant Qn-i- On the other hand, denoting by Yj the j-th column of the matrix Y, 
we obtain from the jump equation ( p.lOD (or, more precisely, from the equation y+ = Y^H^^) 
that 

= -z-"y-(^)( ) 
= -z-"y2_(z) + ,/.(z)z-"y_(z)/(z) 

= -z-"y2-(z)+<^(z)z-"F(z). (2.15) 

Define 

Jiz) = I < 1' 

and consider the 2x2 matrix function 



z-"y2(z), \z\ > 1, 



Z{z) = as (F(z), J{z)) as, ^3 = J ^^j- (2.16) 
The function Z is analytic outside of C, and it has the following asymptotic behavior as z ^ oo: 



For the jump relation on the contour C we have from ( 2.15]) , 



Z+{z) = as (F(z), -yi+(z)) ^3 = ^3 (F(z), z-^Y2^{z) - 0(z)z-"F(z)) .73 



a, {F{z),z~-Y2-{z)) I I ^^f'^ ]<y, = <y, {F{z),z-Y,_{z)) a,a, I J ^^f'^ ) 
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Summarizing the analytic properties of Z we conclude that it solves the following Riemann- 
Hilbert problem: 

• Z is holomorphic for all z ^ C, 

• Z+{z) = Z^{z)S{z), zeC, 
where the jump matrix S{z) is 

= ( ; '""f > ) . (2.18) 

We shall refer to this Riemann-Hilbert problem as Z-RH problem. As in the y-RH problem, 
the solution of the Z-RH problem is unique. Indeed, assuming that Z is another solution, we 
introduce the matrix ratio X := ZZ~^. By the same reasoning as in the case of the y-RH 
problem, we conclude that X is entire. Since 



X{z) = -^I as z 



oo, 



it follows that X = I; and hence, that Z is unique. We note that Y (and hence the resolvent 
R) can be reconstructed from Z using ( |2.16 ). It also should be pointed out that the existence 
of the solution of the Z-RH problem (as well as of the y-RH problem) is equivalent to the 
nondegeneracy of the Toeplitz matrix T„ ((/>), i.e. to the inequality 

Dn{(p) / 0, 

which we always assume. 

Remark. There is a more direct and elegant way to pass to the Z-RH problem which was 
pointed out by the referee of this paper. One first notes that the jump matrix H admits the 
factorization, 

1 j I 1 I [ -I z"" j ' 



H{z) = 

which then suggests the definition 



Y{z) 



Y{z) 'I ^ , \z\< 1, 



(2.19) 



so that the function Y would satisfy the Riemann-Hilbert problem. 
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• y is holomorphic for all z ^ C, 

• ^ I, z^oo, 

*Y_iz) = Y,iz)^l ^7^^)), zee. 

The function Z{z) is related to Y(z) by 

Z{z) = a3Y{z)a3. (2.20) 

We conclude this section by summarizing the relation of the Z-RH problem to the orthogonal 
polynomials on C with respect to the (generally complex) weight 0. This relation is due to 
Deift[^ (see also [|^).0 Let {Pk{z)}k=o i denote the system of the monic polynomials defined 
by 

Pkiz) = z'^ + 
dz 

Pniz)Pm{z)(j){z)— = KKm-, n > 171, 

c 

where bar denotes complex conjugation. Similarly, introduce a second system of polynomials, 
{Pk{z)}k=o,i,...^ by replacing (p with (f) in the definition of P„. Suppose now that 

D;,(0)/O, fc = l,...,n + l. (2.21) 

Then (see |2^) both the sets of polynomials {-Pfe}fc=o,i,...,n and {-P^}fc=o,i,...,n exist, and the 
normalization constants and /i^, k = 1, . . . ,n, are all nonzero. In fact, we have the explicit 
representations 

Dn+l{(t)\z) Dn+l{(l>) ^, , Dn+l{4)\z) ^ Dn+l{^) 

Dn{(P) Dn{(p) Dn{<p) -D„(0) 

(2.22) 

where Dn+i{(j)\z) denotes the Toeplitz determinant Dn+i{4>) whose last row is replaced by the 
row (1 1 z, z , z ). If we define 

Qk = -^PUl/-z)z', (2.23) 

and 

/ P„(z) ^/^P„(/)(/)-0(.')^ \ 
Ziz) =\ \, (2.24) 

V Qn^liz) ^JaQn-lizlizT'^Hz')^, J 

then it is a calculation to show that this Z defines a (unique) solution of the Z-RH problem 
(cf. [|12|, Ql, and [Q). Indeed, the analyticity in C \ C and the jump condition follow from the 
basic properties of Cauchy integral, and the asymptotic condition at z = oo is equivalent to 
the fact that the polynomials P„ and Pn-i are monic orthogonal polynomials with the weights 
(j){z)dz and (j){z)dz, respectively. 

^The Z-RH problem is the analog for polynomials on the circle of the Riemann-Hilbert problem derived in 
for polynomials which are orthogonal with respect to an exponential weight on the line (see also ^). 
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3 Toeplitz Determinants as Integrable Systems 
3.1 Universal recursion relation 

In this section will be an arbitrary continuous function with Fourier coefficients c^j. We 
assume that the associated Toeplitz matrix Tn{4>) is invertible. Then the corresponding matrix 
RH problem is uniquely solvable, and the following equation connects the Toeplitz determinant 
Dn{(t>) with the solution Z of the Riemann-Hilbert problem, 

^ = ^12(0), (3.1) 

where Zjj, i,j = 1, 2, denotes the entries of matrix Z. Indeed, using ( 2.22| ) we have that 

Dn+l 1 



On the other hand, ( |2.24| ) gives 

■^12(0) = —hn, 
zvr 

and ( |3.1| ) follows. 

Remark. One can prove (|3.l| ) using only the connection with the integrable operator Kn{(l>) 
introduced in (p.l|). To see this first note 

- z-z' 2^ 

1 " /zX'^ l-Mz') 



^/ 5Z ( ^/ 



z' ,^ \z' J 27ri 

fc=0 

= Kn{z,z')+h{z)gi{z')^, 

z' 

where / and g are defined in ( |2.3| ) and (p.4|), respectively. Attaching to the functions /, g 
superscript "n" to denote their n dependence, we can rewrite the last equation as an operator 
equation 

Kn+i = Kn + f^®g'^+\ (3.2) 

where the symbol a 6 denotes the integral operator with kernel a{z)b{z'). Recalling the 
definition of F, dJ), it follows from (pi) (cf. pi M) that 



det(l - i^„+i) = det(l - i^„) det(l - [(1 - K^r'm ^ 9i+'] 

= det{l - Kn) det{l - (g) g'l^^) 

= det{l- Kn) {1- trace F^ ^g"^^^) 

= deiil - Kn) (^l - J^F,{z)gi{z)^^ 
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where Fi := and gi := g^. Thus (see also (2.2)) 

Dn+l _ det Kn+l _ , [ ^ r X r 

Dn det Kn 

Recahing (^), we rewrite ( |3.3D as 

Dn+l 



1- / Fi(z)5i(z)-. (3.3) 
Jc z 



which together with (|2.16D yields (|3J1 ). 



>^ii(0), 



3.2 Differentiation formulas 

Here we restrict to the symbol 



AI , 



a=l 



<^(.) = e*^ n ( ) ' (3-4) 



where := —pi^, and we recall (see §1) that ia = max 

M 

a=l 

is the decomposition of the alphabet A into subsets Ai,A2i ■■■■,Am such that pi = pj if and 
only if i and j belong the same Aa- We also recall that 

M 

^ ^ ka = k, 

a=l 

and 

1 >Pi >P2 > ... >Pfc >0, ^Pj = l (3.5) 

j=i 

denote the probabilities assigned to the letters i = 1,2, ... ,k, in the alphabet A. Note that 
from the probabilistic conditions (|3.5|) it follows that 



- 1 < < 0, a = l,...,M, ra^rp,a^P, (3.6) 

and 

5] ko^ra = -1. (3.7) 

The symbols // and (p are related by 

fl{z) = 4>{\lz), 
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and therefore; the corresponding Toephtz matrices are mutually transpose. Thus 

Gi{n;{pi],t) = Dn{4>). 

In what follows, we will write Tn{t), Kn{t) and Dn{t) for Tn{4>), Kn{(j)) and Dn{(j)), respec- 
tively; or Tn{{pi},t), Kn{{pi}, t) and Dn{{pi}, t) if the dependence on pi, . . . ,pk is of interest. 

We shall derive the differential formulas for the Toeplitz determinant Dn{t) with respect to 
the variables t and ra,a = 1,...,M assuming that the latter are subject to restriction (^) 
only, i.e. we only will assume that 

-1 < Tq, < 0, a = 1, M, ra / r^, a / (3. 

The integers kc, will be kept constant. This means that when vary the points ra we do not 



assume restriction (3.1) to hold. We will begin with the t - derivative. 
Since d(j)/dt = zcp, 

at z — z Zttz 

2TTi z — z' 27ri z — z' 27ri 

Let A be the integral operator with the kernel 

Kiz z') = ^''^'^^ - ' ^ 
' z-z' 2m 

Consider the operator product Ai^„: 

{AKn){z,Z) = — / W' 



2ni Jc z — w w — z' 2ni 



2m Jc't^Aw) \z' J ^ ' 2m 2m ^ ' 



27ri 2-Ki 



z ■ 



Recalling the definitions of / and g, (|2.3| ) and( |2.4D , ( |3.8|) can be written compactly as 

^i^n = -/2®52-A(l-i^n). (3.9) 
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From this formula we see that (cf. the derivation of 

— logDnit) = -trace Kl-K„)-i—i^„ 



trace ® 52 + trace A 

F2iz)g2iz)dz, (3.10) 



where we used the fact that 



c 



n I 

trace A = / ctz = 0. 

2m Jc 



RecaUing (^) we convert ( 3.1C| ) into the identity 



d 

— log Dn{t) = -reS;,=oo{Y22{z)), 

which in terms of the Z-function is 

d 

— log Dn{t) = -res^=oo(^;"Z22(z)), 

or equivalently, 

^logZ)„(t) = (ri)22, (3.11) 

where the matrix Fi = Ti{{pi},t) is defined by the expansion, 

OD 



Z{z)=\I + J2^\ z^^^ \z\>l. (3.12) 



■ 1 



Remark. In the basis {z"'}'i^^_^, ( |3^ ) coincides with (3.22) of ^| 



Equation ( 3.11| ) is the t - differentiation formula, i.e. it gives an expression of the t 



derivative of log-D„ in terms of the solution Z of the Riemann-Hilbert problem. We shall 
proceed now with the derivation of the Tq, - differentiation formula. Since (we recall that 
are assumed independent and that the ka are kept constant) 



d ^ ^ kg 
drry z-rr, 



the Tq, - analog of (^]^) reads 



—Kjz,z') = kj-^i^T^ m. 

dva " ' z — z' z' — ra 2m 

kg {z/z'r - 1 1 _ j Jz/zT-l 1 l^iifO _ (3 -^3) 

z — z' z' — rn 27ri 
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Introducing the integral operator with the kernel, 

ka [zjz'Y - 1 1 



we consider again the operator product tv^Kn- The residue type calculations, similar to the 
ones used in the t - case, yield the equation 



(A,K.)(z,z') = |^(l-'A(^')) 



{z/z' 



1 



+ 



{rJz'Y 



1 



Xz' - z){roc - z) {ra - z){ra - z') 
which in turn impies that ( |3.13| ) can be rewritten as 

'r, 



—KJz,z') = 1 - '^(^ ) 

dra " ' liri — z){ra — z') 



+ [K{l-Kn)]{z,z'). 



With the help of the vector functions, 

1 



fiz), ~9{z):-- 



1 



z -rn 



z-ra 

the last equation can be transformed into the following compact form (cf. ( |3.9D ) 

d 



dra 



Kn = kar^f2 <?)9i- kafi <^ 9i + Aa(l - Kn). 



(3.14) 



Let the vector function F{z) = {Fi{z), F2{z))'^ be defined by the equation 

:= (1 - Kj-7„ j = l,2. 

We observe that 



Hz) 



1 



z-rn 



-Y-\Ta)F{z), 



(3.15) 



where the matrix function Y{z) is the solution of the y-RH problem corresponding to Dn{t). 
Indeed by the definition of the vector function F{z) (see ( |2.6D ) its component Fj{z) satisfies 
the integral equation. 



(3.16) 



Fj{z)- / Kn{z,z')Fj{z')dz' = fj{z). 
•J c 

Dividing both sides of this equation by (z — Tq), using the formula 

f{z)9{z') 



Kn{z,z) 



z — z 



and simple algebra, we can rewrite ( |3.16| ) as an equation for the ratio Fj/{ 



z-rn 



Fjiz) 



z-rn 



c 



Kn{z,z') 



z' - rn 



dz! + / f {z)9{z') 
Jc 



F,{z') 



z' - rn 



dz' = fjiz). 
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By applying the operator (1 — Kn) ^ to the both sides of this equation it can be transformed 
into the relation, 



Fjiz) 



z-rn 



+ / F^iz)g{z' 



c 



Z' — Tr 



or 



z 



^F,{z) + Y.F^{z) f m{z')F,iz')- 



dz' = Fj{z), 
dz' 



F,{z). 



The last equation in turn can be viewed as the linear algebraic system for the vector F{z) 

2 



F,iz)-J2A,Aiz) 



1=1 



z-rn 



-Fjiz), j = l,2, 



(3.17) 



where the matrix A is given by the formula, 



A,, = / F,iz')g,{z')- 



dz' 



c 



z' - r„ 



Equation ( 3.15| ) follows directly from ( |3.17| ) in virtue of definition ( p.7| ) of the matrix function 
Y{z). 

We now able to finish the derivation of the Tq, - differentiation formula for the Toeplitz 
determinant In fact from ( ^.14 ) it follows that (cf. the derivation of ( |3.10| )) 



d 

dVa 



log Dn{t) 



( d 
-trace [{l-Kn)-^—Kn 

-fco-r^trace F2®gi + A;atrace Fi® gi — trace Ac 



Krl / F2{z)gi{z)- 



dz 



c 



z-r^ 

where, similar to the t - derivative case, we used the fact that 



+ K / Fi{z)gi{z)- 



dz 



C 



z-rn 



(3.18) 



trace 



nkn 



dz 



0. 



27ri Jc z{z - ra) 

Using now ( ^.15 ) and the fact that dety(2;) = 1 we derive from ( ^.18|) that 

dz 



d 

dra 



log Dn{t) = KrlY2i{r^) f F,{z)gi{z) 

Jc 



KrlYii{rn) / F2{z)gi{z 



dz 



c 



{z - r^y 

+kaY22{ra) I Fi{z)gi{z)—^^^— - kn^Yuira) / ^2(^)51(2) 
Jc [z-raY 



[z - ra) 
dz 



c 



{z - r^y 



or 



_d_ 

drn 



logDnit) 



-fcar2>2i(ra)n'i(^«) + kn,CYii{rn;)Yi^{rn) 
- KY22{ra)Yii{rn,) + knXi2{rn,)Y^^{rn,), 



(3.19) 
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where we use the notation, 



Y'.(r ) — ^^iiMl i j -1 2 



Equation ( 3.19 ) can be also rewritten as 
d 



dr. 



logDnit) = -k^{r2Y2i{ra)+Y22{ra))Yl^{ra) 

+ka{r:^Yn{ra) + Yi2{ra))Y^iK), (3.20) 

which in turn can be transformed into an expression of the dlog Dn{t) / dr^ in terms of the 
function. Indeed recahing formulae ( |2.16| ), (|2.9| ), and (|2.3| ), we see that inside the unit circle 
C the following equation takes place, 

Z(z\-( ^"^11^+^12^ Yn{z) \ 
^'~\-z^Y2i{z)-Y22{z) -Y2i{z)l 



so that (l3^ ) can be converted into the rQ-differentiation formula, 

— \ogDn{t) = -ka[Zn{ra)Z!22{ra) - Z2i{r^)Z[^{r^)). (3.21) 

3.3 Schlesinger equations 

In this section we show that Dn{t) is the Jimbo-Miwa-Ueno r- function of the generalized 
Schlesinger system of nonlinear differential equations describing the isomonodromy deforma- 
tions of the 2x2 matrix linear ODE which has M + 1 simple poles in the finite complex plane 
and one Poincare index 1 irregular singular point at infinity. We will also evaluate the relevant 
monodromy data that single out the Dn{t) from all the other solutions of the Schlesinger sys- 
tem. In the uniform case, when all pi are equal, the system reduces to the particular case of 
Painleve V equation, i.e. we are back to the uniform result of [0O| ]. 

Define 

'^'iz) = eT'^^( I ] (3.22) 



where 



(-z) 



^ /z-r 



a=l ^ 

We note that that the product 
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is our symbol, i.e. the function defined in (^.4|). We also note that is analytic and invertible 
in C \ {0, ri, ...,rM}, and that it satisfies the linear differential equations 

= (3.24) 
'J>?(^) = ^^3<I>°(^), (3.25) 



where is the rational matrix function 
(Subscripts on <I>^ denote differentiation.) Introduce 



$(z) = Z(z)$°(z), (3.28) 

where Z is the solution of the Z-RH problem corresponding to our symbol (p, and consider the 
logarithmic derivative 

B{z) {z)<!>~\z). (3.29) 

The key observation is that B is continuous across the contour C. Indeed, the Z-jump matrix 
S (see ( 2. IS )) admits the following factorization, 

S{z) = $0(z) 



1 1 




} ) (3.30) 



so that the $-jump matrix does not depend on z. In fact we have 

$+(z) = <I>_(z)(^ J \ y zGC. (3.31) 

This implies that 

B+{z)=B^{z), zeC, 

and hence the function B{z) is an analytic function on C \ {0,ri, ...,rA/}. We also recall that 
the only conditions which we impose on the points are the inequalities (p.6|), i.e., 

-l<ro<0, a=l,...,M, r„ 7^ a / /3. (3.32) 

We now calculate the principal part of B at each of its singular points. Since Z is holomor- 
phic and invertible inside of C, it follows from 

B{z) = Z{z)n{z)Z~^{z) + Z^{z)Z-^{z) (3.33) 
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that in a neighborhood of z = r^, 
B{z) ^" 







z — rr 



zir^)i \z-\r^) + j2bjiz-r^y. 



(3.34) 



j=0 



Likewise in a neie; hborhood of z = 0, (|3^ ) and ( |337D imply that 



^(-) = ^^(o)(o J J^-^(o) + E^>^^■ 



(3.35) 



Finally, from ( 3.33| ) and the Laurent expansion ( p. 12 ) we obtain the power series of i? at cxd, 



Biz) = -as + Y.b 



2 — 



Equations (p.34D-(pj.36|) imply that S is a rational function, 



^ t Bo ^ Ba 



with the matrix residues given by 



(3.36) 



(3.37) 



(3.38) 



5„ = -Z(r,)f|j Mz-^r^), a = l,...,M. 



(3.39) 



Thus from ( ^.29 ) we conclude that $ satisfies the linear differential equation, 

$,(z) = 5(z)$(z), (3.40) 

with the coefficient matrix B determined by ( |3.37| )-( |3.39| ). 

Remark. In soliton theory (see [^]), the method that we used to derive (3.40) is called 
Zakharov-Shabat dressing of the vacuum equation ( |3.24|) . We also note that, as is common in 
the analysis of soliton equations, we have moved the exponential factor e^^/^ to the asymptotic 
condition at z = oo. 

Let us now dress the t- vacuum equation ( |3.25| ), i.e. consider the i-logarithmic derivative of 

V{z) := ^t{z)^-\z). (3.41) 
The ^>-jump matrix ( |3.31| ) does not depend on t as well. Hence 



V+{z) = V-{z), zee, 
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and V is analytic on C \ {0, ri, ...,rm}. In fact, since 



V{z) = Z{z)-a^Z-\z) + Zt{z)Z~\z), 



(3.42) 



(cf. ( p.33D ) and Z is holomorphic at the points {0, ri, }, we conclude that V is entire. 
Moreover, from the expansion ( 3.12 ) we have that 



V{z) = -CT3 + -[a^,Ti] + Y,VjZ-^, 



\z\ > 1, 



and hence 



V(z) = ^a3 + ^[as,r,]. 



([L,M] ■.= LM -ML.) 

This in turn yields the t-equation for $, 

$t(z) = V{z)^z), 

where the coefficient matrix V is defined by the equations, 

z 



V{z) 
Vo 



0-3 + Vo, 



k3,ri 



-res. 



iZ(z)z- 



(see also (3.12)). 



(3.43) 

(3.44) 

(3.45) 

(3.46) 
(3.47) 



Equations ( p.40| ) and ( p. 44 ) form an overdetermined system for the function $ in the vari- 
ables z and t. From the compatibility condition, 

we derive the following equation for the coefficient matrices B and V, 

Bt{z)-V,{z) = [V{z),B{z)], (3.48) 

or, taking into account ( p. 45 ), 



Bt{z)-^ = ^[as,B{z)] + [Vo,B{z)]. 



(3.49) 



Since this equation is satisfied identically in z, a comparison of the principal parts of both the 
sides at z = 0,ri, . . . , vm, then leads to the differential relations. 



-Qf = [Y^3 + Vo,B^]. a = 0,l,...,M, 



(3.50) 



(It is notationally convenient to define tq = 0.) 
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The important point now is that the matrix Vq can be expressed in terms of the matrices 
Ba, so that relations (|3.50| ) form a closed system of nonlinear ODEs for the matrix residues 
Ba- In fact, expanding both sides of ( 3.33| ) in a Laurent series at z = oo, using ( 0.27 ), ( 3.37] ), 
and (3.12), and equating the terms of order we have 

m , 

* [rijCJs] + n\ ^ ^ ] +na3 




1 



a=l 




1 



+ na3 




1 



2 

2 iri,cT3j + - era + - I. 



Comparing the last equation with (|3.46) we obtain 



Vn 



a=0 



n n 

0-3 - — /, 



2t 



2t 



so that ( 3.5C| ) becomes 
dt 



2t 



a = 0, 1 



, J., . . . , 



M. 



(3.51) 



(3.52) 



(3.53) 



7=0 



If we vary the points r^, then we obtain M additional linear differential equations for 
$(2) = <P{z,t,ri,...,rM), 



B. 



z-rn 



-$(2), a = l,...,M. 



(3.54) 



Indeed, introducing the - logarithmic derivative, 

Ua{z) := $,jz)cl>-^(z), 

and using exactly the same line of arguments as before, we conclude that Ua is analytic on 
C \ {0, ri, ...,rm}- Simultaneously, the - vacuum equation ( 3.26 ) implies the identity. 



UJz) 



-Z z\ 



Z - Vn 



(0 l^Z-\z) + ZrAz)Z-Hz), (3.55) 



(cf. ( |3.33 ) and ( 3.42| )) which indicates that the only singularity of Ua is a simple pole at 
with 



z = rn 



Z-Tn 



-Z(rn 



Q J ^ Z-\ra) = (see also (^ ) 



as the corresponding principal part. Moreover, taking into account that the asymptotics of 
Z(z) as z ^ 00 does not depend on we conclude that 



Ua{z) — > 0, Z ^ CO, 
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and hence 



Equation ( |3.54 ) now follows. 

The compatibility conditions of equations ( |3.54| ) with ( 3.40 ) lead to the nonlinear r-differential 
equations for the matrices = Ba{t,ri, ...,rM), 



dBg 
dBp 

dB^ 

dra 



[Ba , B^ 



[Bp, Bg] 



a = l,...,M, 



E[Bj , Bee] -I , J 

a = l,...,M. 



(3.56) 
(3.57) 
(3.58) 



which supplement i-equation (|3.53D . 

The total system (( |3.53 ), (3.56)-( 3.58D ) of nonlinear PDEs is the (generalized) system of 
Schlesinger equations which describes the isomonodromy deformations (see e.g. ]l8| , p!9| ]) of the 
coefficients of the 2x2 system of linear ODEs having M + 1 regular singularities at the points 
= 0, ...,M and an irregualr singular point of Poincare index 1 at infinity (see ( |3.40| ), 

(H)), 

d^{z) t ^ B^ 



dz 



B{z)'^{z), B{z) 



t ^ 

a=0 



(3.59) 



The monodromy data of equation ( 3.59| ) which single out the solution of (( p. 531 ), ( p.56D -( p.58D ), 
which we are interested in, coincide, after the proper normalization, with the data of the Z-RH 
problem. More precisely, let us denote ^"^{z) the analytic continuation of ^{z) from \z\ > 1 to 



the whole complex z-plane. Then, the Z-RH problem and equation ( 3.28| ) imply the following 
representations of the function ^°^{z) in the neighborhoods of its singular points. 



$°°(z) 
$°°(z) 
<^^(z) 




, zGUq 



Z G t/oo, ^oo(oo) = /. 



(3.60) 
(3.61) 
(3.62) 



Here ^aiz), ^q{z), and ^od{z) denote the matrix functions which are holomorphic and invert- 
ible in the neighborhoods Ur^, Uq, and Uoo respectively. Formulae (|3.60| )- (|3.62| ) allow us to 
identify the diagonal matrices. 




En 





n + k 



and Er. 




(3.63) 
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as the formal monodromy exponents (cf. |T^, |l9|) of ^°°{z) at the points Tq,, 0, and oo respec- 
tively. The corresponding connection matrices, i.e. the matrices Ca in the representations, 

<^^{z) = ^a{z){z-raf-Ca, a = 0,...,M, 

all are given by 

Ca={l~}], a = 0,...,M. (3.64) 



Since the numbers k^, k, and n are integers, all the monodromy matrices of <l>°°(z) are 
trivial. There are also no Stokes' matrices at the irregular singular point z = oo since the 
asymptotic series ( ^.12 ), as a Laurent series, converges in a disk centered at infinity. Therefore 



the complete monodromy data of the linear system ( p. 59 ) for our random word problem, 



consists of (i) (|3.63| ), the formal monodromy exponents at the singular points, and (ii) ( 3.64 ) 
the corresponding connection matrices. 



3.4 Toeplitz determinant as a r-function 



In this section we shall derive the exact formulae for the logarithmic derivatives of the Toeplitz 
determinant Dn{t,ri, ...,rM) in terms of the matrices Ba which, as we saw in the previous 
section, satisfy the Schlesinger equations (( |3.53| ), ( |3.56 )-( 3.58[ )). To this end we will exploit 
( 3. lip and ( 3.21 ) whose right hand sides we will express via using a technique similar to 
the one that led to (|]52|). We begin with ( |3Tl| ). 

Equation ( |3.52| ) was obtained by expanding both sides of ( 3.33| ) about oo and then equating 
the terms of order z~'^. Let us now analyze the terms of order z~'^. From (3.27) it follows that 



n{z) 



2 



n( 0\_J_^ / 

t ^2 ^ / n 



,-2 



Combining this with expansion 

( CT) of Z, we get the following expression for the order 
term of the right hand side of ( p. 33 ): 

t t n 

- [(73, Ti] Li + - [Fa, as] + [Fi, VLi] + - Ti + - [Fi, ag]. 

The order term of the left hand side follows directly from ( p. 37 ): 

M 

X] '^o.Bo,. 

Q=l 

Equating the two expressions we arrive at 

t t n 

E = 2 [^3,ri]Fi + - [F2,a3] + [Fi,r?i] + - Ti + - [Fi,a3]. (3.65) 
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This equation together with ( |3.51| ) determines Fi in terms of the matrices B^. Indeed, ( 3.51J ) 
gives the off diagonal part of Fi. Using that for L diagonal 

diag[P,L] =0, 

we have from ( 3.65| ) that 

m , 

diagFi = -diag 5^r„B« + -diag([a3,ri]ri) + r!2. (3.66) 

a=l ^ 

Using the identity that for any 2x2 matrix P, 

diag([a3,P]P) = -^[a3,P]V3, 

we obtain from (|3.66| ) and ( |3.51| ) the final expression for the diagonal part of Fi, 

M 1 / Af \ / M \ 

diagTi = -diag ^ r^Ba - - XI ~ - l^^] 51 ~ - -l] (^3 + ^2 

a=l * Vq=0 ^ ^ / \a=0 ^ ^ / 

= -diag raBa - - diag (I XI ~ 2 ~ 2 ) ^ ^"^^ ) ~^ 

M I I n \ 

= -diag X ^o^a - -diag ^ BaB^as - ^"^{(^3 + I)Bacr3 + 1^2 

a=l ya,7=0 a=0 j 

M I ( ^ n \ 

= -diag X ^a^a - - diag X ^aS-,c73 - - X ^"(^3 + /) + f^2. (3.67) 

0=1 ya,7=0 a=0 y 

We also made use of the identities, 

diag (PL) = if diagP = and diagL = L, 

and 

/ M \ 

diag (^E ^« - 2 - 2 = 0. (3-68) 

(The latter follows from ( 3.51| ).) 

We are at last ready to evaluate dlog Dn/dt in terms of B^. To this end it is convenient to 
use 

trace Ti = (which follows from det Z = 1) 
to rewrite ( p. 11 ) in the form 

^logZ)„(t) = -^trace (Ticxs) , (3.69) 
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and then use (l3l67| ) and (|3^ to obtain 



d 1 1 1 ^ 



2t ^ ' 2t 2 

a=0 a, 7=0 a=l 



For future comparison with the Jimbo-Miwa-Ueno r-function it is convenient to use ( p. 681 ) one 
more time and rewrite ( p.70 ) as 

Q ^ M MM ^ M 



0=0 



7=0 



a=l 



j^i=l,2 a=0 

Let us now perform the similar transformations with the right hand side of equation ( 3.21| ). 
We first notice that its subscripts - free form can be written down as 

d 



dra 



log Dn{t) = trace Z-^(ra)Z'(r«)^„ 



(3.72) 



where 





-kn 



is the formal monodromy exponent at (see ( 3.63 )) and in transforming (|]2l|) into (1^721) we 
took into account that detZ(2;) = 1. Secondly, by rewriting equation ( 3.331 ) as the equation 

Z-\z)Z'{z) = Z'\z)B{z)Z{z) - n{z), 

we get the following representation of the product Z~^{ra)Z' {ra)Ea, 

Z-\r^)Z'{r^)E^ = lz'\r^)asZ{r^)E^ + Z-\r^)B,Z[r^)E^ 



2 

t 



7=0 



M 



Ea-J2 + [Z'\r^)Z'{r^)E^,E^]. (3.73) 



(For notational convenience we set, as before, ro := and /cq := —n — k.) 

Using equation (p. 73 ) in the right hand side of equation ( 3.72| ) and taking into account that 

S„ = Z(r„)^„Z-i(r„), 

we arrive to the following r^-analog of ( |3.7[1| ) 



i-Jo,D^it) = ^trace(i?„a3)+E*""^''"''^^ 



7#ai 



M 



kfy t ^ > k' 



■7 



2 Ta 

7#a 



(3.74) 
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Combining equations ( 3.71| ) and ( |3.74| ) we obtain the main result of this section which is the 
fohowing equation for the total differential of the function log Dn{t,ri, ...,rM), 



M 



M 



+ - ^ tr:ace{Baa3)d{rat) 



a,'y—0 



a=0 



^ M M 

+ 9 H (^^a)ij{^B^)j 



dt 

T 



l<i,j<2 0=0 7=0 

-— ^ ^ k^di^rcet) — — ^ ^ k(^k. 



7 



"7 



(3.75) 



a=l 



a, 7—0 



Equation ( 3.75| ) describes the Toeplitz determinant Dn{t) in terms of the solution of the 
Schlesinger system (( 3.53|) , ( ^.56 )-( ^^ )) up to a multiplicative constant (depending on n 
and ka). Simulteneously, this equation shows, upon comparison with the expression (5.17) in 



|18| for the logarithmic derivative of the r-function, that 

M 



TjMU) 



(3.76) 



a ,7—0 



where we use the notation Tjmu = TjMv{t,rQ,ri, ...,rM) for the Jimbo-Miwa-Ueno r-function 
corresponding to the linear system ( p. 591 ) and evaluated for the monodromy data given in ( 3.63| ) 
and dpI). 

Remark. It follows from (|]76D that 

''"jMu vanishes as Tq, — > r-y for some pair (0,7). This 
fact, of course, can be established directly from the definition of the r-function. 



4 Summary of the results 

Recall that Dn{(f>) denotes the Toeplitz determinant associated with the symbol 

(/.(z) = e*^ n (^^) " ' -l<^a<0, a = l,...,M, ra^ri3,a^P, 

and that the generating function G j (n; {pi} , t) is given by the formula 

Gi{n;{pi},t) = Dn{cl)), ra = -pi^. 

We also denote by Dn+i{(p\z) the Toeplitz determinant Dn+i{(p) whose last row is replaced by 
the row (1, z^, z"), and we shall assume that Dn{4>) 7^ 0. 
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The following theorem identifies Dn{<p) as an object of the theory of integrable systems; 
more specifically, as an object of the theory of generalized Schlesinger equations developed in 
11,11]. 

Theorem 1. Let Z denote the 2x2 matrix function defined by 



Z{z) 



/ Dn+i{<t>\z) _J_ r Dn+i{<t>\z') ( j\-n !( J\ dz' \ 



j_ r 

277 JC D„{(l>) 



(4.1) 



where C is the unit circle \z\ = 1 oriented counterclockwise. Introduce the 2x2 matrices 
Ba '■= Ba{t) := Ba{{ra},t), a = 0,1, M , by the equations, 



S„ = -Z(r,)(° ^ ]z-\ra), a = 0,...,M, 



(4.2) 



where 

ro := 0, and ko := —n — k. 
(The invertibility of Z follows from statement 4 below.) Then the following statements hold: 



1. 



M 



M 



d\og Dn{t,ri, ...,rM) 



\ V trace {B^B^) — — ^ + ^ trace{Ba(T3)d{rat) 



Tq, r^ 



a=0 



a ,7—0 

1 dt 

l<i,j<2 a=0 7=0 

-— ^ ^ kctd{rat) — — ^ ] kak'y - 

Tq. r-y 



(4.3) 



a, 7 — 



2. The matrices Ba satisfy the system of nonlinear PDEs (generalized Schlesinger equa- 
tions), 



dBg 
dt 

dBg 
drj 
dBo 
dra 
dBa 
dr„ 



M 



n-tr^ . [B^,Ba] 



2t 

[Ba, B^ 



t 



[Ba,CJ3] + J2 
7=0 

■, a/7 = l,...,M, 



a = 0,l,...,M. 



[Bp, Bg] 

ro-ra' 

[-^7) Ba] 



a = l,...,M. 



(4.4) 

(4.5) 
(4.6) 
(4.7) 
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3. Equations (4-4)^(4i2) compatibility conditions for the system of linear equations 

M 

(Tq -I- 

2 



at 

d^{z) 
drry 



t ^ B, 
■^3 + 2^ 



a=0 

n 



z — r, 
n 



a 



2^^ 2t^^ 2t 



1 \ 



z-rn 



-<I>(z), a = 1 



,...,M, 



(4.8) 

(4.9) 
(4.10) 



which in turn implies that the system (4-4)~(4-''' ) describes the isomonodromy deforma- 
tions of the z- equation (4-^)- 

4- The function Z is alternatively defined as an unique solution of the matrix Riemann- 
Hilbert problem, 

• Z is holomorphic for all z ^ C , 

• Z(z)z-"'^3 ^ /, z ^ oo, 
' 1 z-'^Mz) \ 



Z+{z) = Z^{z) 







1 



zeC. 



(In particular, we have that det Z = 1.) Equation (4_^) can be rewritten in terms of Z as 

M 



dlogDn = 
Also, 

5. The function 



res. 



z'^Z22{z))^dt - ka[Zii{ra)Z'22{ra) " ^'21 (ra)Z{2 (r^)) (ir« . 

a=l 

(4.11) 

Dn+l 



^12(0). 



(4.12) 



M 



z - rr, 



satisfies the linear system (4-S)~(4-lC ) with the matrices Be given by {4-i-)- 

6. The matrices Ba are alternatively defined as the solution of the inverse monodromy prob- 
lem for the linear equation jj.^ characterized by the following monodromy data: 

• the formal monodromy exponents at the singular points r^, oo are given by the 
equations 




-K 



E 



-n 
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• the corresponding connection matrices are 

Ca ~ ^ Q 1 ^ ' ~ ^ 

• the Stokes matrices at the irregular singular point, z = oo, are trivial. 

7. 

Dn{ct>) = e-'^^^'--^- n K-r^\'^ rjuu. (4.13) 



a, '7=0 



where Tjmu = TjMuit,fo,i^i,---,i^M) denotes the Jimho-Miwa-Ueno r-function correspond- 
ing to the linear system ([j-Si ) and evaluated for the monodromy data indicated. Equation 



( 4-15 ) in turn implies the following representation for the generating function Gj{n; {pi},t), 

M 



Gi{n;{pi},t) = e2 [[ \pi^ - Pi^,\ ~ c TjMu{t,0, -pi^, -pij,j), := 0. (4.14) 

Remark 1 . In the uniform case, i.e. when M = 1 and ki = k, the hnear system ( |4.g| ) reduces 
to the 2x2 system of hnear ODEs which has two regular singular points and one irregular 
point of Poincare index 1. In this case, as it is shown in [p!8| , |l^, the isomonodromy equations 
([4.4|) - ([4.7D reduce to the special case of the fifth Painleve equation. Consequently this suggests 
that the uniform generating function Gj{n;t) can be expressed in terms of a solution of the 



fifth Painleve equation. That this is so was obtained earlier in [30| via a direct analysis of the 
Toeplitz determinant Z)„(t). 

Remark 2. The methods developed in this paper can be easily generalized to any symbol 

(j){z) := (f)[z, t) such that 

dz log (/), dt log 

are rational in z. 

Remark 3. In virtue of the Fredholm determinant formula ( p.2| ) for the Toeplitz determinant 
Dn{(p), equation ( ^.13| ) can be interpreted as an example of the general relation |2^ between 
the Jimbo-Miwa-Ueno isomonodromy r-function and the Sato-Segal- Wilson r-function defined 
via an appropriate determinant bundle (see also ||l3| for another example of this relation). 



Remark 4- The generalized Schlesinger system ( [4.4D ~(4.7) appeared earlier in [^] in con- 
nection with the sine kernel Fredholm determinant considered on a union of intervals. The 
corresponding monodromy data, and hence the solution, are different from the ones related 
to the Toeplitz determinant Dn{t). For instance, the sine kernel monodromy matrices are not 
trivial (see [|l^; see also [|l4| for higher matrix dimensional generalizations); in fact, each of 
them equals the identity matrix plus a one dimensional projection. 

Remark 5. From the point of view of the asymptotic analysis of the Toeplitz determinant, 
the most important statement of Theorem 1 is Statement 4. It allows one to apply the Riemann- 



Hilbert asymptotic methods of |11, 
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Remark 6. This paper has been primarily concerned with the isomondromy/Riemann- 
Hilbert aspect of our integrable system. Presumably an analysis of the additional compatibility 
conditions, which arise if one extends (|4.8| )-( pITo| ) by a relevant n-difference equation, would 
lead to a Toda like system, see Okounkov |^5| and Adler and van Moerbeke |3^ . 
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